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Abstract. We analyse available experimental data on the total and differential charged-current cross sec- 
tions for quasielastic Vfj.N and V^N scattering, obtained with a variety of nuclear targets in the accelerator 
experiments at ANL, BNL, FNAL, CERN, and IHEP, dating from the end of sixties to the present day. 
The data are used to adjust the poorly known value of the axial-vector mass of the nucleon. 
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1 Introduction 

A precise knowledge of the cross sections for charged- 
current induced quasielastic scattering (QES) of neutrinos 
and antineutrinos on nuclear targets is a pressing demand 
of the current and planning next generation experiments 
with accelerator and atmospheric neutrino beams, aiming 
at the further exploration of neutrino oscillations, probing 
nonstandard neutrino interactions, searches for proton de- 
cay, and related phenomena. 

The quasielastic cross sections are very sensitive to the 
poorly known shape of the weak axial-vector form factor 
Fa{Q^) of the nucleon. Adopting the conventional dipole 
approximation, this form factor is determined by the axial- 
vector coupling gA = Fa{0) and the phenomenological 
parameter Ma, the so-called axial- vector (dipole) mass 
related to the root-mean-square axial radius by 



/ 2\ 6 
'■a = 

gA 



dFAiQ^) 
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The experimental values of Ma extracted from neutrino 
and antineutrino scattering data and from the more in- 
volved and vastly model-dependent analyses of charged 
pion electroproduction off protons, show very wide spread, 
from roughly 0.7 to 1.2 GeV with tlie formal weighted av- 
erages [TJ[2] 



Ma = 



1.026 ± 0.021 GeV from f^j, experiments, 
1.069 ± 0.016 GeV from tt electroproduction. 



The first value, the common default in most current neu- 
trino simulations, is defined largely by i/^d bubble cham- 
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ber experiments; in many of these experiments, the ex- 
tractions of Ma were based on the naive dipole approxi- 
mation for the vector form factors of the nucleon, along 
with other conjectures. The second value should be in fact 
decreased by about 5%, in order to account for hadronic 
loop corrections (see, e.g., Ref. [1]). 

The results of several selected Vf^d, F^H, and tt^ elec- 
troproduction experiments have been recently reanalyzed 
by Bodek et al. [3], using a new improved description of 
the vector form factors ("BBBA(07)" paramctrization) . 
The obtained world average axial mass is 

Ma = 1.014 ± 0.014 GcV (BBBA(07)). 

This value seems to be in conflict with the new results of 
high-statistics neutrino experiments K2K SciFi [4] (oxy- 
gen target) and MiniBooNE [5] (carbon target), reported 
unexpectedly large while mutually consistent values of the 
axial mass: 



Ma = 



1.20 ±0.12 GeV (K2K), 

1.23 ± 0.20 GeV (MiniBooNE) 



A preliminary analysis of antineutrino data in MiniBooNE 
yields a consistent value of Ma [6]. 

Both K2K and MiniBooNE extractions utilize the up- 
dated vector form factors, from Refs. [7l[8] and [9], re- 
spectively. Within the low-Q^ regions explored in K2K 
and MiniBooNE experiments, the difference between these 
parametrizations and BBBA(07) is comparatively small. 

It can be noted that nuclear effects in the K2K analy- 
sis were accounted within the standard relativistic Fermi 
gas (REG) model [10], while the MiniBooNE analysis used 
REG modified by including an "instrumental" free param- 
eter K which changes the strength of Pauli-blocking. A fit 
of the shape above 0.25 GeV^ (where the variations of 
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K, has no significant impact) leads to an even larger value 
oiMA = 1.25 ±0.12 GeV. 

In this study, which is in a sense complementary to 
that by Bodek et al. [3], we attempt to extract the ax- 
ial mass value by a global statistical analysis of all avail- 
able consistent data on the total and differential QES 
cross sections measured in accelerator experiments with 
and 77,, bcamfQ from ANL rnifnirafTl[^[Tl[T71D^ , BNL 
[19l[20, 21,22,23,24, 25, 26, 27, 28, 29, 30 , FNAL [31,32,33, 
[3T,'35,36,37 , CERN 38 . 39 . 40 , 4 L 42 . 43H44l[45llT6(i47ti48ti 4%! 
[50,51, 52, 53, 54, 55, 56,57,58,59,60,61], and IHEP [Hill 
[64lf65 ,66,67,68,69, 70i(7l| . The detector media used in these 
experiments are hydrogen, deuterium, carbon, aluminium, 
argon, iron/steel, propane, freon, and also propane-frcon 
and neon-hydrogen mixtures. 

In the likelihood analysis, we use the most accurate 
phenomenological parametrizations for the vector form 
factors of the nucleon [72(1 73 j. we take into account all 
known sources of uncertainties, in particular, the system- 
atic errors in the energy spectra of and j/,, beams. For 
description of nuclear effects we apply the standard RFG 
model. We examine possible difference between the values 
of Ma extracted from and data, and cross-check 
our results with the data on distributions measured in 
several experiments. 



2 Quasielastic neutrino scattering off free 
nucleon 

2.1 Structure functions and cross section 

Let us first summarize the well-known phenomenology for 
describing the hypercharge conserved quasielastic reac- 
tions on free nucleon targets 



Here Gf is the Fermi coupling, q ~ k — k' is the four- 
momentum transferred from the incoming (anti)ncutrino 
to the nucleon, = ~q^, Mw is the mass of intermediate 
VF-boson; E^, Ei, P/ ~ ^E\ ~ m'j, and 9i are, respec- 
tively, the incident (anti)neutrino energy, outgoing lepton 
energy, momentum, and scattering angle in the lab frame, 
mg is the lepton mass. The leptonic tensor defined by the 
product of the weak leptonic currents, is given by 

La(3ik, k') = 2 [k'^k/3 + kak'p - gapikk') =F iSal3'y5k''k'^] , 

(3) 

where the upper (lower) sign is for i^i (ve). Assuming the 
isotopic invariance, the hadronic tensor is defined by the 
six structure functions Wi{Q^): 



W^pip,q)=-g^pW^ + ^W2 



2M2 
2M2 



M2 

* 2M2 



We, (4) 



where M is the mass of the "isoscalar" nucleon. Then 
combining Eqs. ^ and (j4]) yields 



al3 



AME,, 



Pf cos ( 



± 



+ 



M 
El, + E( 
M 

Ef, + Pi cos 9. 



Ef- 



Pe COS ( 



2M 



■W2 



2M 



- 2x" 



W3 

(5) 



Vi{k) 



-n{p) 
■P{p) 



■ r{k') 

£+{k') - 



p{p'), 
n{p'). 



(1) 



Here k, k' , p, and p' denote the four-momenta and £ stays 
for e, ^, or r. In this paper, we will neglect the proton- 
neutron mass difference!! since the resulting correction, in 
the Vfj./i'^ case, exclusively works near the reaction thresh- 
old and practically negligible for the energies of our cur- 
rent interest. The general formulas which take this effect 
into account, were derived in Ref. [74] (assuming T and C 
invariance) and in Refs. |751l76j (avoiding these assump- 
tions). 

The double differential cross-section for these processes 
is a convolution of spin-averaged leptonic and hadronic 
tensors Lap and Wap' 



dcJhi. 



GlPi 



dEfd cos Of 



7r(l + QVM^Y 



4ME„ 



(2) 



The Ve 



f-r, and Ut beams from past and current accel- 



erator experiments are not appropriate for measuring the QES 
cross sections. 

^ While our computer code operates with the most general 
formulas and relevant kinematics. 



where >c = mi/2M. 

In order to connect the structure functions with the 
nucleon form factors, we define the charged hadronic cur- 
rent for the QES process (see, e.g., Ref. [77]): 



{P{p')\ja\n{p)) = VudUp{p')ra{p,q)Un{p). 



(6) 



Here Vud is the ud transition clement from the Cabibbo- 
Kobayashi-Maskawa quark-mixing matrix and 



raiP, q) ^laFv + icTaP^FM + Fs 



T-, . Pa I Pa Tji I 9q T-n 1 



(7) 



The form factors Fi are in general complex functions of 
Q^. After standard calculations one finds 



W^{Q^) = 2M^\Vudfw,iQ^)S (2{pq) 



(8) 
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with 

wi =\Fa\^ + x' {\Fv + Fm? + \Fa?) , 

W2 =|Fy|2 + \Fa\'' + X' (|Fa/P + 4|f^Tp) , 

W3 =-2Re[F;(Fy + Fa/)], 

Wi =i [x' {\Fm - 2Fs\ +A\Fp + FtI'') - \Fm\^] 

+ |FsP + ^Re [F* {2Fs - F^) - 2FX{Fp + Ft)] , 

W5 + 2Rc [F^ {Fv - x'Fm) - F^ {Fa - 2x'Fp)] , 
W6 -2Im [F^ (Fv - x'Fm) + F^ (Fa - 2x'Fp)] , 

and x' — /AAI^. The only difference between this result 
and that from Ref. [77] is in the relative sign of the terms 
in which docs not contribute to the QES cross sectionlf] 
Inserting Eqs. dSJ and ^ into Eq. (O gives the com- 
monly known formula for the differential cross section for 
reactions ([T]) on free nuclcon targets: 



'8^(l + QVAf2,)2i?2 



A- 



M2 



B- 



A/2 



C 



A/4 



where 



[251E61E7] (Q2 < 1.2 GcV^) and in the IHEP-ITEP spark 
chamber experiment at Serpukhov [SS] {Q^ < 2.4 GeV^), 
adopting the ad hoc dipolc parameterizations 

Fs{Q^)=isFviO){l + QVMl)-\ 

Ft (0') = ^tFa{0) (1 + Q'A4)~' ■ 

The strongest (but yet not too telling) 90% C.L. upper 
limit for the axial SCC strength has been obtained at 
the BNL AGS experiment as a function of the 
"tensor mass" Mt, assuming conservation of vector cur- 
rent (CVC) (that is = 0), and simple dipole form for 
the vector and axial form factors with My = 0.84 GcV 
and Ma = 1-09 GeV. The limit ranges between 0.78 at 
Mt = 0.5 GcV to about 0.11 at Mt = 1.5 GeV. In so 
much as the contribution of the scalar form factor into the 
QES cross section is suppressed by (m^/A/)^ w 0.01, the 
90% C.L. constraint to the vector SCC strength is even 
less impressive: ^5 < 1.9, assuming = 0, A/5 = 1 GeV, 
and the same My and Ma as above. 

Below, keeping in mind this vagueness, we will assume 
the time and charge invariance of the hadronic current. 
Under this standard assumption, all the form factors are 
real functions of and 

Fq = Ft = 0. 



A =2x'\Fv + FmP - (1 + x')\Fv\'' - x'{l + x')\Fm\ 
+ {l + x')\FA\^-Ax'il + x')\FT\' 

{\Fy^FM\^^\FA^2Fp\^ 

-4(l + x')(|F^P + |Fpp)], 
S = T4x'Rc[Fl(Fy +/^m)] 

+ 4x2Re \F^ {Fa - x'Fp) - F^ {Fy - x'Fm)] , 

C + x'\Fm\'' + \Fa\' + ^x'\Ft\') , 

s={k+p)^ ^2ME^ + M'^, 

u = {k' - pf = - 2MEi, = mj- 2ME^ + Q^. 



2.2 Induced scalar and tensor form factors 



2.3 Vector form factors 

The Dirac and Pauli form factors Fv^m arc related to the 
Sachs electric and magnetic form factors Ge,m'- 



Fv 



Ge + x'G 



M 



1 + x' 



Fm 



G 



M 



Gi 



1 



Isotopic symmetry provides simple relation between Ge,ai 
and elastic electric and magnetic form factors of proton 



and neutron G^" and G^": 



G 



Gp r^n 
M ~ 



\jE — (-Tp; — yj I 



At low Q^, a reasonable description of the electric and 
magnetic form factors is given by the dipole approxima- 
tion: 



The quoted formulas take into account the nonstandard 
G parity violating axial and vector second-class currents 
(SCC) which induce the nonzero scalar and tensor form 
factors Fs and Ft- The most robust restrictions on the 
SCC couplings Fs,t(0) come from the studies of /3 decay of 
complex nuclei (see, e.g., Refs. [751175] and quoted therein 
references). However, these studies are almost insensitive 
to the SCC effects at nonzero . The latter were inves- 
tigated in several (anti)neutrino experiments at BNL 



^ According to Llewellyn Smith, the functions tJs = ojs — l>J2 
and LiJn are, respectively, the real and imaginary parts of a 
unique function. Our examination does not confirm this prop- 
erty for the general case of nonvanishing second-class current 
induced form factors Fs and Ft- 



^E 



Gd, 



j-lpGo, G^E ~ 0; G^j ~ fJ-nGo, 



where Gj 



(1 



A/v 



0.84 GeV, and 



{fin) is the anomalous magnetic moment of the proton 
(neutron). Analyses of the almost all earlier neutrino ex- 
periments were based on this approximation. In this study, 
we utilize two more sophisticated models for the form fac- 
tors GI" and G^f - BBBA(07) [71] and GKex(05)[73]. 

The BBBA(07) model is an accurate Kelly type para- 
mctrization of the current experimental data on the form 
factors G|, G^^^, G^, G^, and ratio G|/G^^, which uses 
the Nachtmann scaling variable 



2 1 
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to relate elastic and inelastic form factors, and imposes 
quark-hadron duality asymptotic constraints at high mo- 
mentum transfers where the quark structure dominates. 
The parametrization is based on the same datasets as 
were used by Kelly [501, updated to include some recent 
experimental results. Quark-hadron duality implies that 
the squared ratio of neutron and proton magnetic form 
factors should be the same as the ratio of the correspond- 
ing inelastic structure functions F2 and F2 in the limit 



2.4 Axial-vector and induced pseudoscalar form factors 

For the axial and pseudoscalar form factors we use the 
conventional parametrizations |77j 



FAiQ^)^FAiO) 1 



2M^ 



Ml 



mi 



(9) 

(10) 



Gn 
M 

r<'p 



Fl 



1 + 4(d/u) 
4 + {d/u) ' 



Here d and u are the partonic density functions. The au- 
thors fit the data under the two assumptions: d/u = and 
d/u = 0.2. One more duality-motivated constraint is the 
equality 

(G"i;/G',V)^ ^ {Gl/Gl,f 

applied for the highest data points for the neutron 
electric form factor included into the BBBA(07) fit. 

The GKcx(05) model is in fact a modification of the 
QCD inspired vector dominance model (VDM) by Gari 
and Kriiempclmann (GK) |81j extended and fine-tuned 
by Lomon [82[|83| in order to match the current and con- 
sistent earlier experimental data. The data set used by 
Lomon includes the polarization transfer measurements, 
which are directly related to the ratios of electric to mag- 
netic form factors, and differential cross section measure- 
ments of the magnetic form factors. The electric form fac- 
tors derived from the Rosenbluth separation of the differ- 
ential cross section are only used for the lower range of 

where the magnetic contributions are less dominant. 
Among several versions of the parametrization considered 
by Lomon, we chose the latest one "GKex(05)" described 
in Ref. [73]. This version incorporates the data that has 
become available since the publication [83| and has a bit 
better x^- The fitted parameters agree with the known 
constraints and the model is consistent with VDM at low 
Q^, while approaching perturbative QCD behavior at high 
. The quark-hadron duality constraint is not imposed. 

Figure [1] shows a comparison of the GKex(05) and 
BBBA(07) parametrizations for the form factors G^" and 
G^^" divided by the standard dipole Gd, against the ex- 
perimental data extracted using either the Rosenbluth 
separation or polarization transfer techniques (including 
a series of double-polarization measurements of neutron 
knock-out from a polarized or ^He targets). The data 
assemblage is borrowed from Refs. [Mi.85..86.,87] and re- 
cent reviews [88ll89] . It is seen from the figure that the 
models are numerically close to each other at low momen- 
tum transfers covered by experiment, but diverge at high 
Q"^. The most serious disagreement between the models 
is in the neutron electric form factor ai ~>2 GeV^. In 
section^ we examine how the model differences affect the 
extracted value of the axial mass. 



where Fa{Q) ~ Qa is the axial coupling, is the charged 
pion mass, and Ma is the axial-vector mass treated as a 
free parameter. In fact, Eq. (|1G|) is a conjecture inspired by 
the hypothesis of partial conservation of the axial current 
(PC AC), expectation that the form factor Fp is domi- 
nated by the pion pole near = 0, and the "technical" 
condition 



1 dFAiQ^) 



FaW dQ^ 



Q2=o 



2< 
All 



< 1, 



which is obviously fulfilled for the experimental lower limit 
of Ma- Since the pseudoscalar contribution enters into 
the cross sections multiplied by {mi/M)'^, the uncertainty 
caused by this approximation may only be important for 
Vt/vt induced reactions (especially in the low-Q^ range, 
see, e.g., Refs. [SnUS]) and it is insignificant for reactions 
induced by electron and muon (anti)neutrinos. 



2.5 Constants 

The most precise determination of Vud comes from su- 
perallowed nuclear beta decays (0"^ 0+ transitions). 
We adopt the weighted average of the nine best measured 
superallowed decays l/j^^^ = 0.97377 ± 0.00027 recom- 
mended by the Particle Data Group (PDG) [12 . Note that 
this value is consistent with that of the PIBETA experi- 



ment at PSI 



V, 



(PIBETA) 



ud 



0.9728 ± 0.0030, obtained 



from the measured branching ratio for pion beta decay 

For the axial- vector and Fermi coupling constants, we 
use the standard PDG averaged values: gA = —1.2695 ± 
0.0029 and Gp = 1.16637 x 10"^ GeV^ [92]. In several 
papers (see, e.g., Ref. [94] and references therein) it is sug- 
gested to use the value G'p = 1.1803 x 10"^ GeV^ obtained 
from 0+ 0+ nuclear /? decays, rather than the standard 
Gp obtained from muon (3 decay. The coupling constant 
G'p subsumes the bulk of the inner radiative corrections. 
However, some neutrino experiments already take the ra- 
diative corrections into account (sometimes in quite differ- 
ent ways) in the measured cross sections. That is why, in 
this study, we simply add the corresponding difference (of 
about 2%) to the overall uncertainty of the fit. Note that 
using the G'p instead of Gf would lead to a few percent 
decrease of the output value of Ma- 
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3 Relativistic Fermi gas model 

Since the main part of the experimental data on the QES 
cross sections for nuclear targets was not corrected for 
nuclear effects, we must take these into account in our 
calculations. In the present work, we use the RFG model 
by Smith and Moniz [TD] incorporated as a standard tool 
into essentially all neutrino event generators employed in 
accelerator and astroparticle neutrino experiments. 

According to RFG, the hadronic tensor Wa/s given by 
Eq. ([4]) must be replaced with the tensor Tafj, which de- 
scribes the bound nucleon. This tensor is of the same 
Lorentz structure as Wa/3 and is defined by the six in- 
variant nuclear structure functions Ti{Q^). Thus, in the 
in the lab. frame 



Tap (Plab,g) 



gapTi + goagopT2 
„<5 



2M, 



H ^TT J 5 



2Mt 

.goaq/B - go go/3 



2Mt 

dp.f{p,q)Wai3ip, q), 



(11) 



where piab = (A^t, 0), Mt is the mass of the target nucleus, 
and 

f{p, q) = v;^l ni{p) [1 - nf{p + q)] . 

The function Hi (p) is the Fermi momentum distribution of 
the target nucleons. satisfying the normalization condition 



J ni{p)dp = 1. 



The factor 1 ~ nf{p + q) (the unoccupation probability) 
takes into account the Pauli blocking for the outgoing nu- 
cleon. The relative velocity v^ci which represents the flux 
of incident particles, is given by 

Vrei = \{kp)\/{E,Mt). 

Explicitly defining the three-momenta q, p, and p, 

9= (0,0,|q|), 

p = (sin6'fe,0,cos6'fe) |q|, 

p = (sin 6p cos 4>p, sin 6'p sin0p, cos 6'p) |p|, 

one obtains 

Wrci = [Ep - \p\ {cos 9k cos 9p + sinOk sin^p sini^p)] /Mt, 
where 



Ep = + 



is the total energy of the bound nucleon and et is the 
effective binding energy. The angle 6k is defined by 



cos Uk 



El+q^ + ml 
2E,\q\ 



For determining the angle ^p, one can use the energy con- 
servation law defined by delta-function 



5{Ep - Ep+q + i^) = 

where v = E„ — Et and 



TH-^jT—fS (cos 6'p - cos 6*°) 



-q2 + 2|p||q|cos6ip + M2. 



is the total energy of the outgoing nucleon. Then the con- 
dition 



cos Or. 



{v + Epf 



{Ep + Ch) 



2N|g| 



must be obeyed. 

The nuclear structure functions are the linear combi- 
nation of the Wi and can be straightforwardly calculated 
from Eqs. (g]) and (fTT|) : 

1 



Ti ^aiWi + -{a2-a3)W2, 



T2 

T3 



99 9 
q - V , ^ , 

^^(a2-a3) + ^a3 



2v 

\q\ 



Wo 



Ml 

Mi 
'\q\ 



M 

Jq\ 



— (3a3 - as) W2 + aiWi + o^^eW^s 
1^ (02 - Sas) + 2a5 



W2 



0.7 



M \ \q\ 
The coefficients are given by 

ai 



as 



f{p, q)dp, 
02 y" f{p,q)p^dp, 

as =^ J f{p, q)p^ cos^ 6pdp, 
a4 =^ J I{p, q)Eldp, 

f{p,q)Ep\p\ cos Opdp, 
.f{p, q)\p\ cos 9pdp, 

a? =^ J f{p, q)Epdp. 

Finally, in order to describe the neutrino scattering off 
a bound nucleon, one should substitute M \ — > Mt and 
Wi I — > Ti in Eq. ([5]); then the differential cross-section 
can be calculated according to Eq. ([2]) (see Ref. [61] for 
more details). Table [T] collects the values of proton and 
neutron Fermi momenta p^" and binding energies e^'" for 
several nuclei, used in our numerical calculations. 



1 

1 
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Table 1. Proton and neutron Fermi momenta and binding 
energies (in MeV) for selected nuclei. 



Nucleus 






Pf 






ZZL 


O K 
ZO.D 


OO 1 


or a 




223 


26.2 


223 


26.1 




225 


26.6 


225 


26.6 


'If 


233 


28.4 


233 


28.3 


?SNe 


230 


27.8 


230 


27.8 




239 


29.5 


239 


29.4 


isAr 


242 


30.7 


259 


35.0 


iFe 


251 


33.0 


263 


36.1 


iBr 


245 


31.5 


270 


38.1 



4 Statistical analysis of the data 

4.1 Description of experimental data 

We have examined and classified all available experimental 
data on quasielastic scattering with AY = 0. Published re- 
sults from the relevant experiments with and beams 
from accelerators at ANL ,11, 12, 13, 14,'T5 l [T6 l fT7 l [T8]. BNL 
[Tni[P,l^,l^ll31[21[inil^,27,28,29,30 , FNAL [3T1,'3?,'^ 
[3a[35,36,37l, CERN r38,39,40,41,42,43l[4ll|45l|46,47,48,4i 
[50':5L52,53,54','55','56,57,,58,59,60,61j, and IHEP 62,63] 
r64.,65..66..67.,68..69.,70.i71j arc included dating from the end 
of sixties to the present day, covering a variety of nuclear 
targets, with energies ranging from about 150 MeV (ANL 
experiments) to about 350 GcV (NuTcV). Pertinent addi- 
tional information was borrowed from the review articles 
and data compilations [Ml[Ml [Wl[MraiOOyiOiyiO'2yi03[ 
I104lll()5lll()6lll()7lll()8lll()9j . 

All the fits are done with the CERN function mini- 
mization and error analysis package "MINUIT" (version 
94.1) [llOj . taking care of getting an accurate error ma- 
trix. The errors of the output parameters quoted below 
correspond to the usual one-standard-deviation (Icr) er- 
rors (MINUIT default). 

For the analysis, we have selected the most statistically 
reliable measurements of the total and differential cross 
sections for each nuclear target, which were not super- 
seded or reconsidered (due to increased statistics, revised 
normalization, etc.) in the posterior reports of the same 
experimental groups. Finally, we include into the global 
fit the data on the total cross sections from Rcfs. 17, 21] 
[22,32,36,37,40,50,55,60,61,611171] and the data for the 
differential cross sections from Refs. [50 l [57 l [64 l l65 p 68 ( [7n 
llOlj . The remaining data are either obsolete, or exhibit 
uncontrollable systematic errors and/or fall well outside 
the most probable range determined through the fit of 
the full dataset; the value of evaluated for each subset 
of the rejected data usually exceeds (3 — 4) NDF. 

Since the differential cross sections da/dQ^ were mea- 
sured, as a rule, within rather wide ranges of the energy 
spectra of and beams, we use only the data from 
such experiments, in which the spectra were known (mea- 
sured or calculated and then calibrated) with reasonably 



good accuracy. All the energy spectra (borrowed from 
Refs. [5Uir5MS51 [Tgnil04mil[[TT^ ) necessary for numerical 
averaging of the calculated differential cross sections and 
distributions were parametrized. To avoid the loss of ac- 
curacy, the precision of these parametrizations was cho- 
sen to be at least an order of magnitude better than the 
experimental accuracy of the spectra themselves. For a 
verification, we have estimated the mean energies of the 
beams for different energy intervals, and have compared 
these against the published values. 

The analyses were performed for neutrino and antineu- 
trino data separately, and for the full set of the v and V 
data together. For each fit, we have included the data for 
either total or differential cross sections, as well as for the 
cross sections of both types together. The main results of 
the analysis are summarised in Tables [5] and [3] and illus- 
trated in Figs.[2HT9l Let us discuss these results in details. 

4.2 Main results of the global fit 

As is seen from Table [21 the differences between the values 
of Ma extracted from the fits of each type, performed with 
the BBBA(07) and GKex(05) models for the vector form 
factors vary between 0.3% and 1.3% that is less than or 
of the order of one standard deviation in the Ma extrac- 
tions and is comparable with the accuracy of the most 
precise measurements of the electric and magnetic form 
factors. The values of x^/NDF are essentially the same 
for BBBA(07) and GKex(05). The differences in the Ma 
values obtained with the two versions of the BBBA(07) 
model corresponding to d/u = and 0.2 (the latter is not 
shown in the table) are less than 0.2% that is practically 
negligible. Therefore, in the following we will solely discuss 
the d/u ~ case. 

The Ma values obtained from the fits to the differential 
cross sections are systematically lower those obtained from 
the total cross sections. The differences amount ~ 1.5% 
(~ 5.7%) for (i^fi) that is (especially in antineutrino 
case) above the statistical error of the fit and is caused 
mainly by uncertainties in the energy spectra of j/^ and 
77^1 and, in lesser extent, in the nuclear effects. 

Figures [2] and [3] show a compilation of the available 
data on the total QES cross sections for the following nu- 
clear targets: hydrogen gT], deuterium pai51ll61fT7[[^ 
|32l|57], carbon [61], aluminium pl[65l[67l[68] . argon [60], 
iron [37], steel [12], propane [43], freon f40U45 | l48 l l50 l l62 ] 
[TOIITTIITUI] . and also propane-freon [52,55, 55] and neon- 
hydrogen ;33,34,3g mixtures. The recent MiniBooNE 2007 
datapoint (carbon target) estimated from the reported 
value of Ma is also shown in Fig. [2] for comparison. 

The compilation does not include obviously obsolete 
data (e.g., ANL 1972 [13], CERN HLBC 1965/1966 [Ml 
[39] ) , as well as the data identical to those reported in the 
posterior publications of the same experimental groups 
(e.g., FNAL 1982 [31], GGM 1978 [S3], IHEP-ITEP 1983 
[66], IHEP SKAT 1986 [69]). The early results of the NO- 
MAD experiment reported in Refs. [551155] . have been con- 
siderably revised (mainly due to corrections in nuclear 
Monte Carlo) [5T] ; the datapoints shown in Fig. [51 are still 
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Table 2. Values of Ma (given in GeV), extracted by fitting the v^^,, 77^, and v^^ +77^ data on total and differential QES cross 
sections, using the BBBA(07) and GKex(05) models for the vector form factors of the nucleon. The x^/NDF values for each fit 
are shown in parentheses. 



BBBA(07) 




GKex(05) 




Ml Ml M"+~ 


Ml 


Ml 




Fit to the total cross sections: 








0.994 ±0.017 1.047 ±0.025 1.011 ±0.014 


0.986 ±0.017 


1.035 ± 0.025 


1.001 ±0.014 


(83/82) (134/62) (220/145) 


(83/82) 


(137/62) 


(222/145) 


Fit to the differential cross sections: 








0.979 ±0.020 0.991 ±0.029 0.983 ± 0.017 


0.976 ± 0.020 


0.982 ± 0.030 


0.978 ±0.017 


(45/48) (26/37) (71/86) 


(45/48) 


(25/37) 


(70/86) 


Fit to the total and differential cross sections: 








0.988 ±0.013 1.023 ±0.018 0.999 ±0.011 


0.981 ±0.013 


1.012 ± 0.019 


0.991 ±0.011 


(128/131) (163/100) (293/232) 


(128/131) 


(163/100) 


(293/232) 



preliminary and are reproduced here by permission of the 
NOMAD Collaboration. 

All the deuterium data quoted in Fig.[2]and freon data 
in Fig. [3] were converted to a free nucleon target by the 
experimenters]! The BNL 1981 experiment [22] had re- 
ported the E,j and dependencies of Ma extracted from 
a fit of the experimental distribution rather than the 
cross section; we quote the BNL 1981 cross section recal- 
culated from Ma by Kitagaki et al. |32]. Similarly, the 
FNAL 1984 rectangle [331JM] and FNAL 1987 datapoint 
|36| were calculated by the experimenters (for free proton 
target) using the Ma value extracted from the measured 

distribution oiv^ events recorded in the Fermilab 15' 
bubble chamber filled with a heavy neon-hydrogen mix- 
ture. The data from several freon experiments (e.g., [40i 
I451I97| ) reported in the original papers in units cm^ per 
nucleon of freon nucleus, were converted to the standard 
units. 

All solid curves shown in the figures were calculated 
using the BBBA(07) model for vector form factors with 
d/u = and always correspond to the best fit value 

Af^+^ = 0.999 ± 0.011 GeV (xVNDF « 1.3) (12) 

obtained from the global fit of neutrino and antineutrino 
data on the total and differential cross sections (see Table 
12]) . We do not show the cross sections calculated with the 
GKex(05) model since the difference will be practically 
invisible. 

The dashed curves in Fig. [2] are calculated with the 
Ma values extracted from the best fit to the (preliminary) 
NOMAD total cross section data alone pT) : 

M'i = 1.05 ± 0.023tat ± 0.07syst GeV, 

Ml = 1.06 ± O.OTstat ± 0.123y.t GeV, ^^^^ 

* The nuclear corrections applied to the deuterium data un- 
der consideration, were treated according to Singh |113j . The 
nuclear effects for the freon data were modeled using a Fermi 
gas approach. 



both agree with the global fit value IT^ . Note that these 
results were obtained with the GKex(05) vector form fac- 
tors. Fitting the NOMAD data with the BBBA(07) form 
factors increases Af^ and A/^ by about 0.8 and 0.9%, re- 
spectively, that still remains well within the errors quoted 
in dni). 

As is seen from the figures, the obtained result, despite 
the non-optimal a-i^d large spread of the data, is not in 
conflict with the main part of the data excluded from the 
global fit. Moreover, it well agrees with the world averaged 
value of 

Ma = 1.014 ± 0.014 GeV, (14) 

obtained in Ref. [3] as a result of their reanalysis of the 
"raw" data from u^d and T7^i? experiments ANL 1973 
dH, ANL 1977 [m, ANL 1982 [H], BNL 1980 [H], BNL 
1981 [22], BNL 1983 [m], BNL 1990 [28], FNAL 1983 
[22], CERN BEBC 1990 [57], and from pion electropro- 
duction experiments after corrections for hadronic effects. 
Note that the values of Ma re-extracted in Ref. [3] from 
each Vfj^d experiment separately spread between 0.97±0.05 
and 1.04 ± 0.06 GeV. It exceeds the difference between 
the results of our analysis of data on total and differen- 
tial cross sections. Both analyses use the same BBBA(07) 
model and mutually supplement each other, since they 
practically do not overlap in the adopted data sets. For- 
mal averaging of the values (|12p and yields 

Ma = 1.006 ±0.009 GeV. 



4.3 Are M'J^ and M^ really different? 

According to the global fit (see Table [2]) , the difference 
between the values of and Af^ obtained by fitting the 
neutrino and antineutrino data separately, reaches about 
3.5% for BBBA(07) and about 3.2% for GKex(05) that 
is above the statistical error in determination of Af^ and 
M^. However, taking into account the systematic differ- 
ence between the fits of total and differential cross section 
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data, as well as high values of x^/NDF, this difference can- 
not be considered statistically significant. Furthermore, 
the fit to the antineutrino data is not stable relative to 
including/excluding some data subsets. In particular, as 
is seen from Fig. [51 the total NuTeV cross sections per 
nucleon bound in iron, averaged over the energy range 
E^,T, = 30 ^ 300 GeV 

a{iy^,n ^-p) = (0.94 ± O.OSstat ± O.OT^yst) x IQ-''^ cm^ 
aiV^p^fi+n) = (l.I2±0.04stat±0.I0syst)xI0~-''«cm2 

(shown in Fig. [2] by rectangles) notably exceed the cor- 
responding best fit curves whereby the NuTeV data [37] 
strongly affects the global fit values of and M^. 

To clarify this point further, we have performed addi- 
tional fits, in which the datasets obtained in experiments 
with non-active targets have been removed. Namely, we 
excluded the highest energy NuTeV total cross section 
data (iron target) [37] and the data on differential cross 
sections measured with the IHEP-ITEP spark chamber 
detector with aluminium filters [551[^51[S5] , since these ex- 
periments do not have an active target to measure recoil 
hadrons and surely remove resonance background. In or- 
der to minimize possible uncertainties in nuclear correc- 
tions, the lowest-energy CERN 1967 total cross section 
data (freon target) [40] were also excluded from these fits. 
The results of this analysis are summarized in Table [31 It 
is seen that the additional reduction of the dataset essen- 
tially decreases the resulting values of Ma- Concurrently 
it improves the statistical quality of the fits to the total 
cross section data, while slightly increases the x^/NDF for 
the fit to the differential cross sections. Besides that, the 
Ma values extracted from the total and differential cross 
sections become bit more consistent. The differences be- 
tween and [-65 MeV for BBBA(07) and -75 MeV 
for GKex(05)] become opposite in sign to those obtained 
from our "default" fit performed with the full dataset. 
However, both M^ and M^ values are still compatible, 
within the la deviation, with the average value of M'^'^ . 
So we may reckon that 

(i) the axial mass extraction is rather responsive to the 
choice of the data subsets and 

(ii) the current experimental data cannot definitely con- 
firm or disconfirm possible difference between the axial 
masses extracted from experiments with neutrino and 
antineutrino beams. 

Similar fit performed for the differential cross section data 
only, from which all the i/^d data were excluded, leads to 
an increase of by about 4.2% (4.4%) for BBBA(07) 
(GKex(05)). However, the statistical error of this fit in- 
creases too. Including into this fit the non-deutcrium data 
on total cross sections diminish the increase of Mj^ to 
about 1.2% for both BBBA(07) and GKex(05). Hence, 
the above conclusions remain essentially unchanged. 

4.4 Further details on differential cross section data 

As is known from the comparison with the low-energy 
electron-nucleus scattering data, the RFG description of 



the low-Q^ region is not enough accurate especially at en- 
ergies below ~ 2 GeV (for recent discussion, see, e.g., Refs. 
|115[|116j and references therein). Moreover, the shape of 
da/dQ^ at < 0.1 GeV is slowly sensitive to varia- 
tions of Ma (see below). Thus, in order to minimize pos- 
sible uncertainties due to nuclear effects, the points with 
< 0.15 GeV^ were rejected from the fit of the differen- 
tial cross section dataset. Leaving these points in the fit 
would lead to a decrease of the output values of M^, M^, 
and A/^'" obtained from the da/dQ^ dataset by, respec- 
tively, 1.8, 3.3, and 2.2% for BBBA(07) and 2.0, 4.0, and 
2.6% for GKex(05) form factors. The corresponding de- 
crease of A/a derived from the full dataset (tr and da/dQ"^) 
is clearly less essential: respectively, 0.7, 1.3, and 0.9% for 
BBBA(07) and 0.7, 1.5, and 1.0% for GKex(05). 

Of course, the mentioned uncertainty still remains in 
the RFG calculations of the total cross sections, since the 
contribution from the low-Q^ region is essential at low en- 
ergies. To illustrate this, we show in Fig. [3] the relative 
contribution of the region < Ql into the total cross 
section, R{Qf) = a (Q^ < Qf) /cr, as a function of Qf, 
evaluated for and V^, QE interactions with carbon at 
several (anti)neutrino energies using Ma = 1 GbYE 

It is 

seen that for neutrino-nucleus interactions R < 0.25 as 
Ql < 0.15 GeV^ and > 0.7 GeV that is for aU energies 
of our current interest. As a result, a few percent error ex- 
pected in da/dQ^ due to inaccuracy of the RFG model for 
the low-Q^ region, becomes nearly negligible in the total 
cross section. However it is not the case for antineutrino 
interactions, for which the ratio R (Qf = 0.15 GeV^) be- 
comes reasonably small {R < 0.3) only for > 2 GeV. 
Therefore the lower energy antineutrino total cross sec- 
tion data may bias an uncontrolled (while still small) ad- 
ditional uncertainty. Fortunately, the major part of the 
data participated in the global fit satisfies the above con- 
ditions and our examination demonstrates that the related 
uncertainty is not weighty. 

Figures [SHZKa) and [SHTTl represent the spectrum-ave- 
raged differential cross sections for several nuclear targets: 
deuterium (Fig. [5]) [57 , aluminium (Figs.[6]and[7|(a)) [63l 
[651111], freon (Figs. [l'[9] and [TO]) [SDUTOllTTl lTOTllTUi] . and 
propane- freon mixture (Fig. [TT]) [551152]. In Fig. [71(b) we 
show (for illustrative purposes only) the axial-vector form 
factor extracted in the IHEP-ITEP spark chamber exper- 
iment [68] • All the quoted data, except those from Ref. 
[70] (superseded by the data from the more recent publi- 
cation by the SKAT Collaboration [71]), model-dependent 
IHEP-ITEP data on Fa{Q^) [68], and a few rejected low- 

datapoints, participate in the global fit. We show the 
cross sections calculated with Ma obtained by individ- 
ual fits to the data of each experiment alone and com- 
pare these against the cross sections evaluated with the 
global-fit value of Ma- All the details are recounted in 
the captions and legends of the figures. The comparison 
demonstrates that the individual and global fits generally 



^ Here a (Q^ < Qi) is defined as an integral of da/dQ^ from 
the kinematical minimum of to ~ Qi- 
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Table 3. The same as in Table [2] but after exclusion of the datasets from experiments with non-active targets (NuTeV 1984 
[37], IHEP-ITEP 1981,82,85 [MMIISI]) and the lowest-energy data of CERN 1967 gO] (see text for details). 



BBBA(07) 




GKex(05) 




Ml Ml M^+~ 


Ml 


Ml 


M'f 


Fit to the total cross sections: 

0.986 ±0.021 0.855 ±0.046 0.958 ± 0.019 

(42/52) (38/35) (88/88) 


0.977 ±0.021 

(42/52) 


0.837 ± 0.046 
(38/35) 


0.948 ±0.019 

(89/88) 


Fit to the differential cross sections: 








0.966 ±0.024 0.971 ±0.042 0.967 ± 0.021 

(33/33) (16/22) (49/56) 


0.963 ± 0.024 

(34/33) 


0.959 ± 0.043 
(15/22) 


0.962 ± 0.021 

(49/56) 


Fit to the total and differential cross sections: 

0.977 ±0.016 0.912 ±0.030 0.962 ± 0.014 

(75/86) (58/58) (137/145) 


0.971 ±0.016 

(76/86) 


0.896 ± 0.031 
(57/58) 


0.954 ± 0.014 
(138/145) 



do not contradict each other. The differences are within 
the experimental errors and arc not of systematic nature. 

As a further test of the global fit, we show in Fig. 
[12] the flux-weighted differential cross sections da{v^n 
fj,~p)/dy and daiv^p — > fi~^n)/dy (divided by energy), 
which were measured with the Gargamelle bubble cham- 
ber filled with liquid freon and exposed to the wide-band 
CERN-PS i^^ and beams. Several analyses of these data 
samples are available from the literature (see Refs. [46] 
ITTlliMTUT] and also Ref. [TOS] for a review). Figure [H] 
shows two representative versions taken from Refs. |47j 
and |101j - the preliminary and final results of the GGM 
experiment, respectively. The data are shown for the five 
narrow instrumental ranges: 1 — 2, 2 — 3, 3 — 5, 5 — 11, 
and 5 — 20 GeV. The measured cross sections were con- 
verted from freon to a free nucleon target by the experi- 
menters, after accounting for Fermi motion of the nucleons 
and Pauli suppression of quasielastic events. 

For a qualitative comparison, we have performed indi- 
vidual fits to the GGM data, separately for neutrino and 
antineutrino differential cross sections. In order to reduce 
possible error introduced by RFG calculations of nuclear 
effects, the energy range of 1 — 2 GeV has been excluded 
from this likelihood analysis. As is seen from the figure, 
the Ma value extracted from the neutrino subsample does 
not contradict to that from the global fit, while it is not 
so for the antineutrino data subsample where the discrep- 
ancy is essential. This discrepancy can be attributed (at 
least, partially) to the vagueness of the model for nuclear 
effects used in the analyses of the GGM data. Since the 
details of the GGM nuclear Monte Carlo are not available, 
we do not include this data sample into the global fit. We 
note, however, that the inclusion of these data (also with- 
out the low-energy datapoints) into the fit only leads to 
a small decrease of the output values of M^, M^, and 
Mf - by, respectively, 0.4, 2.2, and 0.9% for BBBA(07) 
and 0.3, 2.0, and 0.8% for GKex(05) form factors. The 
corresponding x^/NDF values remain nearly the same. 



4.5 distributions 

An additional fruitful set of available data is the dis- 
tributions dN/ dQ^ of the QES events measured in several 
experiments with different nuclear targets. Usually just 
dN/dQ^ is considered as the observable most appropriate 
for extracting axial mass value, since it is less dependent 
of the flux and spectrum uncertainties in comparison with 
the differential or total cross sections. However, in compar- 
ison with the differential cross section, the distribution 
has two drawbacks: it contains an uncertainty due to nor- 
malization, and it is generally less responsive to variations 
of Ma at high Q^. Figure[T3l illustrates the second point. It 
shows the distributions and differential cross sections 
for Vfj^ and quasielastic scattering off free nucleons, eval- 
uated with different values of Ma and normalized to the 
corresponding quantities calculated with Ma = 1 GeV. 
The calculations are done with the fixed values of energy 
corresponding to the mean (anti)neutrino beam energies 
in experiments [Mll^51[5Sll62j . It is seen from the figure 
that the region Q2 < 0.15 GeV^ strongly affected by the 
nuclear effects, is sensitive to Ma for dN/dQ^ and less sen- 
sitive for da/dQ^\ the situation is opposite for the high 
region for which the nuclear corrections are less important. 

We use the measured distributions for a consistency 
test of our analysis. For illustration, we show the four 
sets of data on distributions measured in experiments 
HLBC 1969 (propane) [H] (Fig. [H]), IHEP SKAT 1981 
(freon) [Sg (Fig. [H]), CERN GGM 1979 (propane-freon 
mixture) [SS] (Fig. [Tl]), and FNAL E180 (neon-hydrogen 
mixture) [33l[34] (Fig. [17]) . The curves shown in the figures 
are calculated with the global-fit Ma and normalized to 
the data after fitting of the normalization factor N . The 
shaded bands indicate the uncertainty due mainly to in- 
determination of this factor. The obtained best-fit values 
of N should be compared with these evaluated directly 
from the experimental data (all values are shown in the 
legends of the figures). One can see that the agreement is 
excellent everywhere. So, we may conclude that this test 
was quite successful. 
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Another important confirmation of our result is a rea- 
sonably good agreement with the Ma value extracted in 
our earlier analysis of the data on total inelastic v^N and 



liV CC cross sections and relevant observables |117| . 



Finally, Fig. [TH] presents a comparison of the total QES 
cross sections for i/e, i^^, ^^e, I'/j, and Vr interactions 
with free nuclcons, calculated with the obtained best-fit 
value of Ma = 0.999 ± 0.011 GcV by using the BBBA(07) 
model of vector form factors. The shaded bands reproduce 
the uncertainty due to the la error in Ma- 



5 Discussion and conclusions 

We performed a statistical study of the QES total and dif- 
ferential cross section data in order to extract the best-fit 
values of the parameters Ma- Our main results are sum- 
marized in Table [2] are, of course, model dependent and 
can be recommended for use only within the same (or nu- 
merically equivalent) model assumptions as in the present 
analysis. The best-fit values of the axial mass obtained 
by different fits do not contradict to each other and agree 
with the recent re-extraction of Ma from v^d, F^tH, and 
pion clcctroproduction experiments, reported in Ref. 
They are also in agreement with the preliminary result of 
high-statistical NOMAD experiment at CERN, as well as 
with the numerous earlier data which were not included 
into the likelihood analysis. It has been demonstrated that 
removing the data subsets obtained in experiments with 
non-active targets, particularly the NuTeV dataset, leads 
to a further decrease of the extracted values of Ma (see 
Table |3]). In other words, there is no way to increase the 
Ma value which follows from essentially all (anti)neutrino 
data on total and differential QES cross sections. 

On the other hand, our best-fit value of Ma is in a 
conflict with the mean values of Ma reported by K2K and 
MiniBooNE Collaborations IHIS] , even after accounting for 
the maximum possible systematic error of our analysis 
related primarily to its susceptibility to the choice of the 
data subsets. To expound the problem, let us consider the 
representative K2K result with more details. 

The Ma value reported in Ref. [J has been obtained 
with a water target by fitting the distributions of muon 
tracks reconstructed from neutrino-oxygen quasielastic in- 
teractions by using the combined K2K-I and K2K-IIa data 
from the Scintillating Fiber detector (SciFi) in the KEK 
accelerator to Kamioka muon neutrino beam. The exper- 
imental data from the continuation of the K2K-II period 
were not used in the analysis of Ref. j4j . The best-fit values 
of Ma obtained from the K2K-I and K2K-IIa data subsets 
separately are, respectively, 1.12 ± 0.12 GeV (x^/NDF = 
150/127) and 1.25 ± 0.18 GeV (xVNDF = 109/101). 

Figure [12] shows the Vf^n — > ^,~p total cross section 
per neutron bound in oxygen, recalculated from the fitted 
values of Ma derived in Ref. from the distribution 
shape for each reconstructed neutrino energy. It is neces- 
sary to underline here that the authors do not consider 
their result for each energy bin as a measurement^ but 
rather a consistency test. All calculations represented in 
Fig. [19] were done with our default inputs that introduces 



an uncertainty of at most 2%; this uncertainty is added 
quadratically to the quoted error bars. Also shown are the 
cross sections evaluated by using our best fit value p2]). 
the K2K value of 1.20±0.12 GeV, and the value of 1.1 GeV 
used as a default in the recent neutrino oscillation analyses 
to the data from K2K |118[|119j and Super-Kamiokande I 
|120| . A significant systematic discrepancy is clearly seen 
at > 1 GeV. Since the energy region covered by the 
K2K analysis extends to about 4 GeV, it seems problem- 
atic to explain this discrepancy by the inapplicability of 
the REG model alone. 

Considering that the low-energy K2K and MiniBooNE 
data are in agreement with each other and do not contra- 
dict to the high-energy NuTeV results, we may conclude 
that the new generation experiments for studying the qua- 
sielastic neutrino and antineutrino interactions with nucle- 
ons and nuclei are of urgent necessity, in order to resolve 
the inconsistencies between the old and new measurements 
of the axial-vector mass. 
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Note on the recent SciBar result 

In the recent paper by Espinal and Sanchez [121j FI the 
nucleon axial- vector form factor has been determined from 
an analysis of neutrino-Carbon interactions in the K2K 
detector fully active SciBar tracking calorimeter. The best 
fit value of the axial-vector mass obtained in this analysis 
from the distribution of events and using the BBA(04) 
vector form factors [122] is 



M 



(SciBar) 
A 



1.144 ± 0.077 (fit)I^;^;^ (syst) GeV 



(with xVNDF = 17.2/9, 8/9, and 9.8/9 for, respectively, 
1-track, 2-track QE, and 2-track non-QE events). It is in 
agreement with Ma previously measured at SciFi detector 
in K2K [4] . Formally (within the quoted errors) the SciBar 
result does not contradict to the world average value 
but does not support it. 



^ Which became available to us when this paper has already 
been published 
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Fig. 1. Comparison of the GKex(05) and BBBA(07) models for the electric and magnetic form factors of proton and neutron 
(divided by the standard dipole Go) with the data from electron scattering experiments. The data compilation is taken from 
Refs. [84ll85ll86ll87ll88ll89) . The two versions of the BBBA(07) parametrization are shown for the neutron form factors. 
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Belikov et al., IHEP-ITEP 1981 (Al) 
Belikov et al., IHEP-ITEP 1982 (Al) 
Belikov et al., IHEP-ITEP 1985 (Al) 



Fig. 2. Total quasielastic u^n and i/f^p cross sections measured in experiments with deuterium, hydrogen, carbon/propane, 
aluminium, and iron/steel targets at ANL 1969 [12], ANL 1973 [H], ANL 1975 [IS], ANL 1977 |16II17| . BNL 1980 [2l], BNL 
1981 [22], FNAL 1983 [32], FNAL E180 1984 [33ll34] (rectangle in top right panel), FNAL E180 1987 [36], NuTeV 2004 [37] 
(points and rectangles in bottom panels), CERN HLBC 1969 [43], CERN BEBC 1990 [57] (points and rectangle in top left 
panel), CERN NOMAD 2008 [61] (preliminary), IHEP-ITEP 1981 [63], IHEP-ITEP 1982 [65], and IHEP-ITEP 1985 [571168]. 
The deuterium and neon-hydrogen data were converted to a free neutron/proton target by the authors of the experiments. The 
MiniBooNE 2007 point [5] recalculated from the reported value of Ma = 1.23±0.20 GeV is also shown for comparison. The error 
bars represent the total errors which include the flux normalization uncertainties. The solid curves and narrow shaded bands 
are calculated with the BBBA(07) model for the vector form factors, with Ma ~ 0.999 ± 0.011 GeV, the value obtained from 
the global fit to a subset of the full data set of total and differential cross sections (233 data points) . The points shown by grey 
symbols are excluded from the fit, being either superseded by newer experiments, or not satisfying our selection criteria. The 
dashed curves and corresponding bands are the cross sections obtained by fitting the NOMAD 2008 alone with the GKex(2005) 
vector form factors (separately for i/^ and 77^ data). 
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Fig. 3. Total quasielastic //^n and v^j^p cross sections measured with the freon and propane-freon filled bubble chamber ex- 
periments CERN HLBC 1966 [39], CERN HLBC 1967 gO], CERN GGM 1973 [45], CERN GGM 1975 [48l[97], CERN GGM 
1977 [50], CERN GGM 1978 [52], CERN GGM 1979 [MliSS], IHEP SKAT 1981 [12, IHEP SKAT 1988 [TQ, IHEP SKAT 1990 
[71| . and IHEP SKAT 1992 |104) . The point recently obtained in experiment with the Liquid Argon Time Projection Chamber 
(LAr TPC 2007) [60] is also shown. The SKAT datapoints were converted from freon to a free neutron/proton target by the 
authors of the experiments. The error bars represent the total errors which include the uncertainties due to flux normalization 
and nuclear Monte Carlo. The solid curves and narrow shaded bands are calculated with the BBBA(07) model for the vector 
form factors, with the global fit value of Ma = 0.999 ± 0.011 GeV. See caption of Fig. [2] for more details. 
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Fig. 4. The ratio R = a {Q'^ < Ql) /a vs. Qi, evaluated for /^^^ and //^ quasielastic interactions with carbon target at several 
(anti)neutrino energies. The Ma value is taken to be 1 GeV. 
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Fig. 5. Flux-weighted differential cross section for u^n fi^p measured in the WA25 experiment with the CERN bubble 
chamber BEBC filled with deuterium and exposed to high-energy beam at the CERN-SPS [57] ■ The data were converted 
to a free neutron target by the authors of the experiment. The curves are the calculated cross sections averaged over the 
experimental energy spectrum borrowed from Ref. [112) . The energy range and estimated mean energy are given in the 
legend. The dashed curves are for the best fit to the WA25 data, while the solid curves correspond to the global fit to all QES 
data. Shaded band represents la deviation from the best-fitted value of Ma given in the legend. 
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Belikov et al., IHEP-ITEP 1982 
Target: Al 

Event number: 422 + 39 
3 < fi^ < 30 GeV, <£„ > = 7.1 GeV 
=0.999 +0.011 GeV (global fit) 
M, =0.946 +0.070 GeV (best fit) 
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Belikov et al., IHEP-ITEP 1981 
Target: Al 

Event number: 284 + 38 
3 < £^ < 30 GeV, > = 6.2 GeV 
=0.999 ±0.011 GeV (global fit) 
M, =0.905 +0.087 GeV (best fit) 




Fig. 6. Flux-weighted differential cross sections for i^,in —^ fi^p (a) and F^p — > /i+n (b) measured in the IHEP-ITEP experiment 
with a spark chamber detector with aluminium filters and exposed to the U70 broad-band and 77^ beams of the Serpukhov 
PS |63ll65j . The inner and outer bars indicate statistical and total errors, respectively; the overall systematic error of about 
10% is due mainly to uncertainties of the flux normalization and scanning/triggering efflciencies. The curves are the calculated 
cross sections averaged over the experimental Ufj, and energy spectra borrowed from Refs. |68lll04j . The energy range and 
estimated mean energies are given in the legends. The dashed curves are for the best fit to the IHEP-ITEP data, while the solid 
curves correspond to the global fit to all QES data. The points shown by grey symbols are excluded from the fits (see text). 
Shaded bands represent la deviations from the best-fitted values of A'Ia given in the legends. 
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Fig. 7. Flux-weighted semisum of differential cross sections for u^^ii — > fi^p and V^p fi'^n (a) and axial-vector form factor 
Fa{Q^) (b) measured in the IHEP-ITEP experiment with a spark chamber detector with aluminium filters and exposed to the 
U70 broad-band Vfj, and beams of the Serpukhov PS [68) . The error bars represent the total errors which include the overall 
systematic error of about 10% (due mainly to uncertainties of the flux normalization and scanning/triggering efficiencies). The 
curves in panel (a) are the calculated semisum of the cross sections each averaged over the experimental and 77^ energy 
spectra borrowed from Refs. [6811104) . The energy range and estimated mean energy are given in the legend. The dashed curve 
is for the best fit to the quoted IHEP-ITEP data, while the solid curve corresponds to the global flt to all QES data. Shaded 
bands in panels (a) and (b) represent la deviations from the best-fltted values of A'Ia given in the legends. 
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Bonetti et al., CERN GGM 1977 
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=0.999 ±0.011 GeV (global fit) 
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Fig. 8. Flux-weighted differential cross sections for //^n — > /i p (a) and v^j^p - 
chamber Gargamelle filled with heavy freon and exposed to the CERN-PS 



and 



n (b) measured with the heavy-liquid bubble 
'u beams 



[SOIIlOlj . The error bars contain the 
statistical fluctuation and the indetermination on the and fluxes. The curves are the calculated cross sections averaged 
over the experimental and 77^ energy spectra given in Ref. [50]. Only the events with E^jj > 1.5 GeV were accepted. The 
dashed curves are for the best flt to the GGM 1977 data, while the solid curves correspond to the global flt to all QES data. The 
points shown by grey symbols are excluded from the fits (see text). Shaded bands represent la deviations from the best-fitted 
values of Ma given in the legends. 
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Fig. 9. Flux-weighted differential cross sections for Uf^n — > (a) and 77^p — > fj.'^n (b) measured with the freon filled bubble 
chamber SKAT exposed to the U70 broad-band Vfj, and beams of the Serpulihov PS [7011104] (see also Refs. [B^ for the earlier 
analyses of the same data sample). The data were converted to a free nucleon target by the authors of the experiment. The inner 
and outer bars indicate statistical and total errors, respectively; the systematic error includes the uncertainties due to the cross 
section normalization and nuclear Monte Carlo. The curves are the calculated cross sections averaged over the experimental v^i 
and 77 J, energy spectra borrowed from Ref. [104] . The energy range and estimated mean energies are given in the legends. The 
dashed curves are for the best fit to the SKAT 1988 data, while the solid curves correspond to the global fit to all QES data 
(the SKAT 1988 data are excluded from the global flt). Shaded bands represent la deviations from the best-fitted values of Ma 
given in the legends. 
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Fig. 10. Flux-weighted differential cross sections for i/f^n fi^p (a) and F^^p fi'^n (b) measured with the freon filled bubble 
chamber SKAT exposed to the U70 broad-band i/fj, and 77^ beams of the Serpukhov PS [71]. The data were converted to a free 
nucleon target by the authors of the experiment. The inner and outer bars indicate statistical and total errors, respectively; the 
systematic error includes the uncertainties due to the cross section normalization and nuclear Monte Carlo. The curves are the 
calculated cross sections averaged over the experimental i/f^ and energy spectra borrowed from Ref. |104j . The energy range 
and estimated mean energies are given in the legends. The dashed curves are for the best fit to the SKAT 1990 data, while the 
solid curves correspond to the global fit to all QES data. The points shown by grey symbols are excluded from the fits (see 
text). Shaded bands represent la deviations from the best-fitted values of A4a given in the legends. 
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Fig. 11. Flux-weighted differential cross sections for Vfj,n ^Tp (a) and 77^p — > ^'^n (b) measured with the bubble chamber 
Gargamelle filled with light propane-freon mixture and exposed to the CERN-PS i^^ and 77^ beams [55ll52j . The inner and outer 
bars in panel (a) indicate statistical and total errors, respectively; the error bars in panel (b) contain the statistical fiuctuation 
and the indetermination on the flux. The curves are the calculated cross sections averaged over the experimental Vf_i and F^ 
energy spectra given in Refs. [S^ and [SS], respectively. Only the events with E^jj > 1 GeV were accepted. The dashed curves 
are for the best fit to the GGM data, while the solid curves correspond to the global fit to all QES data. The points shown 
by grey symbols are excluded from the fits (see text). Shaded bands represent la deviations from the best-fitted values of Ma 
given in the legends. 



22 



K. S. Kuzmin et al.: Quasielastic axial-vector mass from experiments on neutrino-nucleus scattering. 



2.6 
2.4 
2.2 

2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 

0.6 

0.4 

0.2 


2.6 

2.4 

2.2 

2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 





1 - 2 GeV 




2-3 GeV 




3-5 GeV 



I 5 -11 GeV 

i 




: I 5-20 GeV 



0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 

y y 



0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 1 

V y y 



- Haguenauer et al., CERN GGM 1974 
• Musset and Vialle, CERN GGM 1978 
Target: CFjBr 



= 0.999 ± 0.011 GeV (global fit) 
= 0.933 ± 0.050 GeV (best fit, v^) 
M, = 0.776 ± 0.075 GeV (best fit, v ) 



Fig. 12. Flux-weighted differential cross sections {1/E,^)da{v^n — > fi p)/dy (top panels) and {1/ E^)da{Ufj,p —^ fi'^n)/dy 
(bottom panels) measured with the heavy freon filled bubble chamber Gargamelle exposed to the wide-band CERN-PS and 
beams [iTlfTCTT] . The data from Refs. [13 (range 5-11 GeV) and [lOl] (ranges 1 - 2, 2 - 3, 3 - 5, and 5 - 20 GeV) represent 
two different analyses of the same data sample (see also Refs. |46II49|[T05] for other versions). The measured cross sections were 
converted to a free nucleon target by the authors of the experiment. The quoted error bars are the total errors which include 
the uncertainties in the i/^ and fluxes and nuclear Monte Carlo. The curves are for the calculated cross sections averaged 
(for each energy range indicated in the panels) over the experimental and energy spectra taken from Ref. [lOlj . The 
dashed curves correspond to the Ma values obtained by fitting the GGM 1978 data from the energy ranges 2 — 3, 3 — 5, and 
5 — 20 GeV and GGM 1974 data from the range 5 — 11 GeV (separately for neutrino and antineutrino cross sections). The 
range 1 — 2 GeV is excluded from the analysis in order to minimize the error in modelling the nuclear effects. The solid curves 
correspond to the global flt to all QES data (the GGM data are not included in this fit). Shaded bands represent la deviations 
from the best-fitted values of A'Ia given in the legend. 
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Fig. 13. The distributions dN/dQ^ and differential cross sections dajdQ^ vs. for v^^n and F^p quasielastic scattering, 
calculated with different Ma = 0.85, 0.90, 0.95, 1.05, 1.10, and 1.15 GeV and normalized to the corresponding quantities 
calculated with Ma = 1 GeV at four fixed values of energy corresponding to the mean (anti)neutrino beam energies in 
experiments HLBC 1969 [43], Gargamelle 1979 [56], SKAT 1981 [62], and FNAL 1984 [33ll34] (see Figs. [TiHTTl below) . The 
curves in the four upper panels end up at the kinematical boundaries. 
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Fig. 14. Flux-weighted distribution for Ufj,n fi^p measured with the CERN heavy-liquid bubble chamber (HLBC) filled 
with propane and exposed to the CERN PS i^^ beam [43] • The curve is the distribution calculated with Ma obtained from the 
global fit, averaged over the experimental i/^ energy spectrum from Ref. [lllj . and normalized to the HLBC 1969 data. The 
spectrum is estimated to be accurate within ±15% (the error includes an estimate of systematic effects). The energy range and 
estimated mean energy are given in the legends. Shaded band represents la variation from the average due to uncertainties in 
Ma and normalization factor A''. 
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Fig. 15. Flux-weighted distribution for — > fi^p measured with the freon filled bubble chamber SKAT exposed to 
the U70 broad-band Vfj, beam of the Serpukhov PS [62]. The data were converted to a free nucleon target by the authors of 
the experiment. The inner and outer bars indicate statistical and total errors, respectively; the systematic error includes the 
uncertainties due to the fiux normalization and nuclear Monte Carlo. The curve is the distribution calculated with A/a obtained 
from the global fit, averaged over the experimental energy spectrum from Ref. [104) . and normalized to the SKAT 1981 data. 
The energy range and estimated mean energy are given in the legends. Shaded band represents la variation from the average 
due to uncertainties in Ma and normalization factor A*'. 
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A^= 39.6 + 2.3 (exp) 
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Fig. 16. Flux-weighted distribution for V^p — > fi'^n measured with the bubble chamber Gargamelle filled with light 
propane-freon mixture (87 mole per cent of propane) and exposed to the CERN-PS 17^ beam [56]. The error bars contain both 
statistical and systematic errors. The curve is the distribution calculated with Ma obtained from the global fit, averaged over 
the experimental Vf^ energy spectrum from Ref. [56], and normalized to the GGM 1979 data. The energy range and estimated 
mean energy are given in the legends. Shaded band represents la variation from the average due to uncertainties in Ma and 
normalization factor A*". 
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Fig. 17. Flux-weighted distribution for 77pp — » fi'^n measured in the FNAL E180 experiment with a 15' bubble chamber 
filled with heavy neon-hydrogen mixture (64% of neon atoms) and exposed to the FNAL wide-band 77^ beam [331134) (see also 
Ref. [31] for an earlier version). The curve is the distribution calculated at the mean antineutrino energy of 12.7 ±0.2 GeV, with 
Ma obtained from the global fit and then normalized to the E180 data. [The spectrum averaging procedure cannot be applied 
here, since the spectrum has been evaluated just from the quoted distribution.] Shaded band represents la variation 
from the average due to uncertainties in Ma and normalization factor A''. 
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Fig. 18. Total quasielastic cross sections for electron, muon and r neutrino and antineutrino interactions with free nucleons 
calculated with the best-fit value of Ma = 0.999 ± 0.011 GeV using the BBBA(07) vector form factors. Shaded bands represent 
the uncertainty due to the Icr error in Ma- 
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Fig. 19. Comparison between the QES Vfj, cross sections per neutron bound in oxygen, evaluated with several values of the 
axial mass. The solid curve with narrow band is calculated with our best fit value of Ma; the dashed curve with wide band 
corresponds to the K2K extraction of Ma [J ; the dash-dotted curve is calculated with the current K2K and Super-Kamiokande I 
default Ma = 1.1 GeV [11911120] . The points represent the K2K cross section reconstructed (with our version of RFG model 
and BBBA(07) vector form factors) from the best-fit values of Ma extracted for the five energy bins, as quoted in Fig. 9 of Ref. 
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